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Abstract 

The influence of the boundary layer properties on the heat trans- 
port in turbulent thermal convection is experimentally investigated in 
a cell with a rough bottom plate. It is shown that the standard 2/7 
exponent of the convective heat flow dependence on the Rayleigh num- 
ber, usually observed in cell with smooth boundaries, increases if the 
roughness has power law distributed asperity heights and the thermal 
boundary layer thickness is smaller than the maximum asperity size. 
In contrast a periodic roughness does not influence the heat transport 
law exponent. 

PACS numbers: 47.27.Te, 05.40. +j, 42.25-p 
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In the last years many models and experiments have been done in order 
to understand the heat transport properties of turbulent thermal convection 
in a fluid layer heated from below, that is Rayleigh Benard convection [I], 
H [| |], H 0- These properties are characterized by the dependence on 
the Raleigh number Ra of the non-dimensional heat flow, that is the Nusselt 
number Nu. In many experiments it has been observed that, for Ra > 10 6 , 
Nu has a power law dependence on Ra, that is: Nu = aRa 1 ' . Specifically, 
in fluids with Prandtl number of about 1, a and 7 take the following values 
for 10 6 < Ra < 10 11 : 7 = 2/7 and a ~ 0.2 @, |, |. 

In order to construct a reliable model of the heat transport law several 
aspects of the turbulent flow have to be solved. For example it is still un- 
clear how the heat transport is influenced by the size distribution of thermal 
fluctuations in the boundary layer and by the coupling of these fluctuations 
with the mean circulation flow (MCF). Thermal fluctuations of the boundary 
layer are associated with thermal plumes || [| and the MCF is a large scale 
convective roll involving all the cell containing the convective fluid || |{J. 

The role of MCF on the heat flux has been studied in several experiments, 



but the values of 7 and a are not modified by either the perturbation [10, 11 



or the suppression [12] of the MCF. No influence on the value of 7 has been 



observed in numerical simulation where no-slip boundaries were used [13 
In two recent experiments [13, ITS] boundaries with periodic roughness have 
been used. In this case a is much larger than the value measured in cell 
with smooth boundaries, that is the heat flux is enhanced. However it is 
important to stress that a periodic roughness does not modify the value of 
7 which is still 2/7 as in the case of smooth boundaries. Finally it is worth 
to mention that in a recent experiment the transition toward the ultimate 
regime has been observed Q. This transition, which manifests itself with 
an increasing of 7 for Ra > 10 11 , has been explained by the change of the 
dissipation properties in the thermal and viscous boundary layers (see also 
ref.[0). One of the consequences that one may extract from all these exper- 
iments and simulations is that the value of 7 could be mainly controlled by 
thermal fluctuations (the thermal plumes) and by their size distribution in 
the boundary layer. Therefore a perturbation of this distribution may change 
the value of 7. 

The purpose of this letter is to show that using a very rough boundary, 
with power law distributed asperity heights, it is possible to strongly modify 
the value of 7 = 2/7 observed in cell with smooth boundaries. The influence 
of the boundary roughness on transport properties is a very important topic 
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in turbulence [ITR O. This topic has not been widely studied in turbulent 



thermal convection. In refs. JTj, [15| only a periodic roughness has been con- 
sidered. As already mentioned, this kind of roughness does not modify the 
value of 7 but only that of a. Therefore, in this letter we want to stress the 
difference between a periodic and a random roughness. 

The experimental apparatus has been already described in ref.[^] and we 
recall here only the main features. The cell has horizontal sizes L x = 40cm 
and L y = 10cm and two different heights d equal to 20cm and 10cm. With 
these two cells filled with water (at an average temperature of 45°C corre- 
sponding to a Prandtl number, of about 3) we are able to cover the interval 
10 7 < Ra< 10 10 . The bottom copper plate plate is heated with an electrical 
resistor. The top copper plate is cooled by a water circulation and its tem- 
perature is stabilized by an electronic controller. All the apparatus is inside 
a temperature stabilized box. The temperature of the plates is measured in 
several locations. Local temperature measurements, of the turbulent flow, 
are done with two small thermocouples (PI, P2) of 0.04cm in diameter with 
a response time of 5ms. The probes PI and P2 are located at (L x /4, L y /2) 
and at (L x /2, L y /2) respectively. Both probes can be moved along z with 
micrometric devices in order to measure the mean temperature profile and 
that of the temperature fluctuations as a function of Ra. To measure the 
heat flow we first estimate the heat losses of the cell ||. These heat losses 
are then subtracted from the heating power to evaluate the fraction of heat 
effectively transported by the convective water. 

We perturb the bottom boundary layer by changing the roughness of the 
bottom plate. This roughness is made by small glass spheres of controlled 
diameter glued on the bottom copper plate, with a very thin layer of thermal 
conductive paint. We use N sets of spheres such that each set j is composed 

by spheres having the same diameter hj, with 1 < j < N and hi < h 2 < 

Hn- The sphere number P{hj) in each set is selected in order to produce 
a well defined power law distribution, that is P(hj) = A h~^. Here A is a 
normalization factor such that the ensemble of spheres covers uniformly the 
copper plate surface, that is L x x L y = tt / AJ2j=i,N h]P(hj). The spheres are 
mixed and randomly glued on the bottom plate. The roughness properties 
can be changed by modifying N, £, the minimum sphere diameter hi and the 
maximum sphere diameter h^. 

The roughness has two important characteristics lengths hi and h^. 
These lengths have to be compared with one of the main characteristic length 
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of turbulent thermal convection that is the thermal boundary layer thickness 
A = d/2/Nu. Indeed if A >> hjsr or A < hi it is reasonable to think that 
no effect on the convection thermal properies will be observed. In contrast 
if hi < A < /ijv several important changes could be produced. The mean 
roughness height h does not seem to play any important role. 

In order to understand the role played by the roughness on the boundary 
layer, we performed several experiments. In four experiments, labeled I, II, 
III, IV respectively, the surface roughness had power law distributed asper- 
ities. A periodic roughness, with only one characteristic length, was used 
in another experiment labeled V. Specifically the roughness parameters, in 
the different experiments, took the values indicated in Table I. The results 
of these five experiments have been compared to those obtained in the same 
cell with smooth boundaries || [L2|]. 

The non-dimensional convective heat flow Nu versus Ra measured in a 
cell with smooth boundaries is compared in fig.l with that measured in the 
experiments (III and IV). We clearly see that in the three cases Nu is a power 
law function of Ra, that is Nu cx Ra 1 with 7 ~ 2/7 in the smooth plate case, 
7 ~ 0.45 in experiment IV (£ = 1) and 7 = 0.35 in experiment III(£ = 2). 
It is important to notice that in the experiments III and IV the thermal 
boundary layer thickness A is always smaller than h^ in all the interval 
10 8 < Ra < 10 9 5 where these two experiments have been performed. Indeed 
from fig.l one sees that at Ra = 10 s , Nu ~ 16, therefore A = 0.63cm < h^ in 
experiments III and IV. These measurements clearly show that, when A < h^ 
has a value comparable to that of the roughness, the dependence of Nu as 
a function of Ra is strongly modified and 7 is a function of £, that is the 
exponent of the roughness height distribution P(h). Specifically 7 increases 
when the roughness height distribution becomes more flat. 

In order to show that for A > h N no effect on Nu is observed we describe 
the results of the experiments (I, II). These two experiments have the same £ 
of experiment III but the important difference is that A becomes comparable 
to h n in the middle of the Ra spanning range. Specifically from the Nu 
measured in these experiments, plotted as a function of Ra in figs. 2a) and 
2b), we find A ~ 0.2cm at Ra ~ 1.5 10 9 for experiment I and A = 0.2cm at 
Ra ~ 3 10 8 for experiment II. In fig. 2a) we clearly see that when A < h^, 
that is for Ra > 10 9 the dependence of Nu changes and we find 7 = 0.35 as 
in experiment III. The value of 7 in this figure and in the next is certainly 
not very precise because of the very limited scaling range. However the exact 
value of 7 is not very important for the discussion. What we want to show 
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here is just the clear change of trend for Ra > 10 9 . In contrast for Ra < 10 9 , 
that is for A > h^, we see that the experimental points are parallel to those 
corresponding to the smooth plate. In fig.2b) the results of experiment II are 
directly compared with those of experiment I. Experiment II has the same 
roughness of experiment I and a smaller d therefore A < at Ra ~ 3 10 8 . 
Indeed we see that Nu begins to increase faster for Ra > 3 10 8 . At the 
same time we notice that in both experiments I and II all the points for 
which A > /itv are aligned on the same straight line parallel to that of the 
smooth case. Thus we see that in order to have an interaction of the thermal 
boundary layer with the roughness A should be smaller than h^. 

To show that this is actually the case we have measured the profile of 
the temperature and of the temperature fluctuations as a function of z in 
experiment I for two values of Ra. These profiles are plotted in fig. 3 as a 
function of 0.5,2/A = zNu/d. For A > the profile with roughness is very 
close to that with a smooth plate. In contrast for A < the dependence 
as a function of z of the temperature fluctuation rms and of temperature 
is strongly perturbed by the presence of the roughness, which induces the 
appearance of a second maximum in the rms profile. The position of this 
maximum is close to h^. Thus the profile shape confirms that the dependence 
of Nu as a function of Ra is modified by the roughness only when > A. 

Finally we compare the results of the experiments I,..,IV with those of 
experiment V which has a periodic roughness. The curve Nu versus Ra 
measured in the experiment (V) is plotted in fig.4. The presence of a jump 
in the curve is clearly observed. The location of the jump corresponds to the 
value of Ra where A ~ hw — hi. However above and below this jump the 
slope of the curve is very close to that with smooth plates. These results agree 
with those of ref [|14], where a periodic rough plate was used. They also tell 
us that when A < hi the roughness does not influence the value of 7. From 
the comparison of the results of experiment V with those of experiment III 
and IV one deduces that an important ingredient for modifying the exponent 
7 is the presence of a roughness with a power law distributed height and that 
a periodic roughness does not influence the value of 7. 

In spite of these important changes in the behaviour of Nu versus Ra 
in experiments I, II, III and IV the bulk properties are not modified by the 
presence of the roughness. The histograms and spectra of the local temper- 
ature fluctuations measured in the center of the cell by probe P2 and on the 
side by probe PI are the same with and without roughness. Furthermore 
the frequency f c of the slow oscillation which is related to the MCF period, 
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has the same dependence on Ra observed in experiments with smooth plates 
and Pr ~ 0(1) || [12|, |J]. We find that our data are compatible with a law 
f c ~ x/d 2 Af Ra 0A9 with Af ~ 0.06 without roughness and A/ ~ 0.05 with 
roughness. This result is in agreement with those of ref . || where it is shown 
that f c is always proportional to Ra 0A9 even for Ra > 10 11 where 7 > 2/7 
as in our experiment ref. Q]. Therefore our experiment and that of ref. 0] 
seems to indicate that there is a negligible influence of the value of 7 on the 
dependence of f c on Ra and on the statistical properties. 

At the moment we are unable to construct a model which explains the 
dependence of 7 on £. The prediction of a recently proposed model does not 
agree with our observations |18]| because this model predicts an increasing of 7 
for increasing £ and we observe just the contrary. Nevertheless one can try to 
understand why the presence of a random roughness is so important. From 
refs. |T(]|, [T2 ] we know that any perturbation of the MCF does not change 
the heat transport. Furthermore several models do not need to rely upon the 
MCF in order to explain the Nu cx Ra 2 / 7 law || 0. Therefore if one assumes 
that the heat transport is mainly due to plumes, whose characteristic size is 
close to A, then if A >> it is obvious that no perturbation of the heat 
transport is observed and everything goes as in the smooth case. In contrast 
when h\ < A < and the roughness has a power law distribution, the 
plumes and the boundary layer cannot construct their characteristic length 
because the roughness contains many different lengths. As a consequence 
the modified A changes the dependence of Nu versus Ra. Finally if A < hi 
no influence on 7 is observed because the system is locally equivalent to a 
smooth one. It is also important to recall that in ref. fij the increase of 7 
for Ra > 10 11 has been justified by the change of the dissipation properties 
in the boundary layer. From all of these observations one may extract the 
very important conclusion that the turbulent heat transport is dominated 
by the thermal fluctuations close to the boundary layer and by their size 
distribution. If this distribution is perturbed by the presence of a random 
roughness the transport properties are modified too. Any realistic model 
of thermal convection should take into account these results, which could be 
very useful in the study of convection in geophysical flows where the presence 
of smooth plates is certainly a very idealized case. 

We acknowledge very useful discussion with R. Benzi, D. Lohse, J.F. 
Pinton and E. Villermaux. 
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Experiment 


P{hj) 


N 


h\ [cm] 


h]s; [cm] 


h[cm] 


d[cm] 


I 


h- 2 


3 


0.06 


0.4 


0.08 


20 


II 


h- 2 


3 


0.06 


0.4 


0.08 


10 


III 


h- 2 


5 


0.06 


1 


0.08 


20 


IV 


h- 1 


5 


0.06 


1 


0.08 


20 


V 


5{h — h N ) 


1 


0.2 


0.2 


0.2 


20 



Table 1: Roughness parameters: experiments I, II, III and IV have power law 
distributed asperity heights whereas experiment V has a periodic roughness 
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Figure 1: Dependence of Nu on Ra: with smooth bottom plate (+) and 
with a rough bottom plate in experiment III (*) with £ = 2 and in experiment 
IV (o) with £ = 1. . In experiments III and IV A was always smaller than 
h N . 
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Figure 2: Dependence of Nu on Ra: with smooth bottom plate (+) and 
with a rough bottom plate in experiment I(o) and in experiment II (x). The 
roughness exponent was £ = 2. a) The boundary layer thickness is smaller 
than /ijv at Ra > 10 9 There is a clear change of slope when A < h^. b) 
The results of experiment I are compared with those of experiment II. For 
experiment II, A < Iin for Ra > 3 10 8 which corresponds to the transition 
point . 
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Figure 3: Vertical profiles of the mean temperature (a) and of the rms of the 
temperature fluctuations (b). The symbol (+) corresponds to measurements 
done with a smooth plate whereas (o) and (*) correspond to the rough case 
with \> h N (o) (Ra = 3 10 8 ) and A < h N (*) (Ra = 10 10 ). 
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Figure 4: Dependence of Nu on Ra: with smooth bottom plate (+) and 
with a rough bottom plate in experiment V (o). The roughness has just one 
size h-N = h\ = 2mm in this case. 
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